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Biased transport of elastic cytoskeletal filaments with alternating polarities by 
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We present a simple model for the bidirectional dynamics of actin bundles with alternating polar- 
ities in gliding assays with non-processive myosin motors. In the model, the bundle is represented 
as an elastic chain consisting of monomers with positive and negative polarities. The motion of the 
bundle is induced by the pulling forces of the underlying motors which stochastically attach to the 
monomers and, depending on the polarity of the monomers, pull them in the right or left direction. 
We demonstrate that perfectly a-polar chains consisting of equal numbers of monomers with positive 
and negative polarities may exhibit biased bidirectional motion with non-zero drift. This effect is 
attributed to the elastic tension developed in the chain due to the action of the myosin motors. 
We also show that as a result of this tension, the attachment probability of the motors is greatly 
reduced and becomes strongly dependent on the length of the chain. These surprising effects point 
to the necessity of considering the elasticity of the cytoskeleton in theoretical studies of cooperative 
dynamics of molecular motors. 

PACS numbers: 87.16. Ka, 87.16. Nn, 87.16. Uv, 87.16.A- 
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Motor proteins are molecular machines that convert 
chemical energy into mechanical work by ATP hydroly- 
sis. They "walk" on the microtubule and actin cytoskele- 
ton and pull vesicles or organelles across the cell [1]. The 
intracellular transport of cargoes is achieved mainly by 
the action of individual motors which propagate along 
the cytoskeleton tracks in a direction determined by the 
intrinsic polarity of the filaments [3]- Other processes, 
such as cell motility and mitosis, require the cooperative 
work of many motors. Muscle contraction, for instance, 
involves the simultaneous action of hundreds of myosin 
II motors pulling on attached actin filaments and caus- 
ing them to slide against each other Q. One interest- 
ing outcome of cooperative action of molecular motors 
is their ability to generate bidirectional motion [4] . Bidi- 
rectional movement results from the competition between 
two populations of motors that work against each other 
in opposite directions ,4-8,] . The direction of motion flips 
from one direction to the other due to stochastic events 
of binding and unbinding of motors to the filament which 
tip the force balance between the two motor groups. 

The dynamics of motor-filament systems are often 
studied using in vitro motility assays in which the fila- 
ments glide over a dense bed of immobilized motors and 
their motion is tracked by fluorescent microscopy Q . Re- 
cently, we used such a motility assay to study the dynam- 
ics of actin bundles induced by the cooperative action of 
myosin II motors [lO| . The bundles in these experiments 
were composed of short actin segments which, through 
a sequence of fusion events, assemble into fllaments with 
randomly alternating polarities. Such a-polar bundles ex- 
hibit "back and forth" bidirectional motion. We showed 
that the distribution of "reversal times" (i.e., the dura- 
tions of unidirectional intervals of motion) take an expo- 
nential form P{At) ~ exp (— At/Tj-cv), where Tj-cv is the 
characteristic reversal time of the bidirectional motion 
[lOl . [ll| . Detailed analysis of the dynamics of many bun- 



dles revealed that t^cv is of the order of a few second and 
has no dependence on the length (number of monomers, 
N) of the bundle. This result was in marked contra- 
diction with previous theoretical models which predicted 
that Trov grows exponentially with N [4] . To resolve this 
disagreement, we have introduced a model that takes into 
account the elastic energy stored in the actin bundle due 
to the action of the working motors [I^, UM ■ The elas- 
tic energy modifles the rates at which motors attach to 
and detach from the actin and eliminates the exponential 
dependence of t^cv on N. 

Our previous theoretical treatment of cooperative bidi- 
rectional motion was based on a mean fleld calculation 
of the actin elastic energy, ignoring both: (i) the sequen- 
tial order of the polarities of the monomers, and (ii) the 
positions along the filament where the pulling forces of 
motors are applied. The mean field elastic energy scales 
a,s E ^ NNc, where Nc is the number of attached mo- 
tors [lO|, [il|. Within the mean field picture, the bidi- 
rectional motion on perfectly a-polar tracks consisting of 
an equal number of monomers with right-pointing ( "pos- 
itive") and left-pointing ("negative") polarities has no 
bias, i.e., the intervals of motion in both directions occur 
with equal probability. In this letter we discuss an in- 
teresting effect related to the elasticity of the actin. We 
show that a-polar elastic filaments may exhibit a biased 
bidirectional motion and achieve a net migration along 
the motors-coated surface. For myosin Il-actin systems, 
we find that the drift velocity is typically 2-3 orders of 
magnitude smaller than the velocity of a single myosin 
II motor and is comparable with the speed by which the 
motors move the a-polar bundle cooperatively during the 
intervals of unidirectional motion. This newly identified 
mechanism of propagation may, therefore, be relevant to 
processes of active self-organization of cytoskeletal struc- 
tures during which filaments are transported and joined 
with each other by motor proteins. 
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FIG. 1: A schematic drawing of the system: A chain of con- 
sisting of N monomers with alternating polarities, connected 
to each other by identical springs. The chain lies on a "bed" 
of motors, some of which are connected to the monomers. A 
connected monomer with positive (negative) polarity feels a 
pulling force of size +/ ( — /). Disconnected monomers expe- 
rience no force. 



To demonstrate the effect, we consider the chain illus- 
trated in Fig. [11 consisting of N monomers connected by 
{N — 1) identical springs with a spring constant k. Each 
monomer may be either free and experience no pulling 
force (/ = 0), or attached to one motor in which case it 
is subjected to a force of magnitude / which is directed 
to the right (+/) for monomers with positive polarities 
and left (— /) for monomers with negative polarities. The 
moving velocity of the filament is given by ^ = /totai/'^j 

where /total — J2i=i fi i^ ^^'^ snm. of motor forces ap- 
plied on the monomers and A is the friction coefficient 
of the chain. A chain of TV monomers has 2^ connec- 
tion configurations, where each such configuration can be 
represented by a vector C of size N specifying the state 
(connected/disconnected) of each monomer. For exam- 
ple, a chain of 4 monomers in which the first and third 
monomers are connected to motors will be represented by 
C = (1,0, 1,0). Let us also introduce a vector S whose 
components are related to the polarities of the monomers. 
The vector 5 = (1, 1, —1, 1), for instance, corresponds to 
a chain of 4 monomers in which the polarities of the first, 
second, and fourth monomers is positive while the third 
monomer has a negative polarity. The drift velocity can 
be calculated by averaging over all possible connection 
configurations of the motors (all possible values of the 
vector C): 
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where Pj is the occurrence probability of the configura- 
tion, and the subscript j has been added to A to account 
for possible variations in the friction coefficient between 
the different configurations. The probability Pj depends 
on (i) the number of attached motors in the configura- 
tion, Nc{j) = \\Cj\\^, (ii) the attachment probability of 
a single motor, q, and (iii) the total elastic energy of the 
springs Ef: 
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where /? = {kBT)~^ is the inverse temperature and 
Z is the partition function of the system. The elas- 
tic energy is the sum of the energies of the springs, 
Ef — X]j=i ^?/2fc, where Fi is the force stretching (or 
compressing) the z-th spring. The forces Fi can be cal- 
culated using the following steps: (i) calculate the mean 
force / = ftotax/N = f{C ■ S)/N, (ii) calculate the access 
forces acting on the monomers /* = fCiSi — /, and (iii) 
sum the access forces applied on all the monomers located 
on one side of the spring F, = X^Li f* = " E/Ij+i /*■ 
Our analysis is based on the assumption that variations 
in C (which occur when motors attach to or detach from 
the actin track) lead to instantaneous changes in the ve- 
locity of the filament which should always be propor- 
tional to the total exerted force. This assumption is ex- 
pected to hold for low Reynolds numbers where inertia 
can be neglected. 

Let us analyze the dynamics of a chain of size A^ = 4. 
There are six different a-polar sequences for a chain of 
this length: Sia = -Sib = (1,1,-1,-1), S2a = ~S2b = 

(1, -1, 1, -1), Ssa = -S-ib = (1, -1, -1, 1)- It is easy to 
prove that the drift velocity (see Eqs. ([T]) and ([2])) van- 
ishes for the first four sequences which are antisymmet- 
ric with respect to refiection around the midpoint. This 
is not the case with the last two symmetric sequences. 
To see this, consider the sequence S^a and assume, for 
simplicity, that Xj = A. In the limit q <^ 1, one can 
ignore the configurations in which more than one out of 
the four monomers is connected to a motor. This leaves 
us with only five configurations: (i) Cj — (0,0,0,0), for 
which V-i = and Pj = (1 - g)VZ, (ii) Cj = (1,0,0,0) 
and Cj ^ (0,0,0,1), for which Vj = //A and Pj = 
g(l -g)3e-(7/8)W/2fe)/Z, and (iii) Q = (0,1,0,0) and 
Cj = (0,0,1,0), for which Vj = -//A and P^ = q{l - 
g)3g-(3/8)(,3/V2fe)/^^ Substituting this in Eq. ^ gives 
V^MSsa) - -V^MSsb) ^ -2(//A)g[e-(3/8)(«V2fe) _ 
g-(7/8)(/3/^/2fc)]^ For /3p/2k < 1 we find that the drift 
velocity increases with a third power of the motor force, 
Vdrift — — (/5/2fcA)/'^. This power law has a different ex- 
ponent than 1 - the scaling exponent for the velocity of 
stiff polar chains. 

To further investigate this effect, we calcu- 
lated the drift velocity for chains of A^ = 4M 
monomers with sequences of the form S = 
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2M 



M 



(— 1, . . . , — 1, 1, . . . , 1, — 1, . . . , — 1). Our results are 

summarized in Figs. [2][A) and (B). Fig. [2][A) is based 
on a calculation in which the friction coefficient (see 
Eq. ([T])), Xj — XqN, while in Fig. [2jB), we assumed 
that Xj = XoNciJ). The resuhs for A^ < 28 have 
been derived using a full statistical calculation of the 
partition function, while for larger A^ they have been 
obtained from Monte Carlo simulations. The model 
parameters were assigned the following values which 
are representative of myosin Il-actin systems |lll . uM- 
PP/2k = 0.002, q = 0.1, and //Aq = 6 ^m/sec. Both 
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FIG. 2: The drift velocity Vdrift as a function of the length 
of the chain. The friction coefficient A-, is proportional to the 
number of monomers A'^ in (A) and the number of connected 
motors Nc{j) in (B). The lines are guides to the eye. 



figures show that for small chains of size N < 200, the 
drift velocity increases rapidly with TV. For larger chains 
{N > 200), Vdrift behaves differently in Figs. E^A) and 
(B). In the former it decreases with N, while in the 
latter it saturates and increases again for N > 600. 
Note also that the different scales of the y-axis in both 
figures. These differences can be attributed to the 
different values of Xj used in the cases represented by 
Figs. mjA) and (B). Since for each configuration, the 
ratio between the friction coefficients in both cases 
rx = Xf/>^f = Nc{j)/N < 1, the drift velocity in 
(B) must always be larger than in (A). Fig. [3] depicts 
the mean value of rx (i.e., the mean fraction of con- 
nected monomers) as a function of N . For N < 50, 
(rx) ~ g = 0.1 and, accordingly, the ratio between the 
drift velocities in (B) and (A) in this regime is close 
to one order of magnitude. For N > 50, (rx) drops 
to values much smaller than q, which implies that the 
friction coefficient per monomer decreases with N in 
case (B) and explains why the drift velocity remains high 
and does not decrease sharply as in (A) . The decrease in 
the mean fraction of connected monomers can be traced 
to the fact that configurations with larger Nc{j) have, 
in general, higher elastic energies and, therefore, smaller 




FIG. 3: The mean fraction of monomers connected motors, 
(Nc/N) as a function of A'^. The solid line represents the mean 
field result Eq.© with c = 0.75. 



statistical weights. For the mean field elastic energy, 
Ef/ksT = c{(ip/2k)NN,{j) [lS[lll, one gets 
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where c is a dimensionless constant of the order of 1. 
For c = 0.75 and N < 1000, this expression (solid 
line in Fig. [3]) gives a fair agreement with the com- 
putational results. For larger values of N (i.e., when 
{Nc/N) becomes very small), the expression tends to 
overestimate the rate of decrease in the mean fraction 
of connected monomers (or, equivalently, the effective 
attachment probability). The decrease in the attach- 
ment probability of the motors is another, indirect, man- 
ifestation of cooperativity between the motors which is 
mediated through the forces that they jointly exert on 
the actin track. Eq. ([3]) suggests that the elasticity of 
the track can be neglected for small filaments whose size 
N ^ (/3/2/2fc)^^ = N* . In this regime, the two coopera- 
tivity effects discussed here which are associated with the 
elasticity of the actin filaments disappear: (i) The drift 
velocity 1/drift - (/3/2A:A)/3 = {f/X){N*/N) « (//A) is 
vanishingly smaller than the typical speed by which the 
bidirectionally moving bundle propagates in each direc- 
tion, and (ii) the fraction of attached motors {Nc/N) ~ q 
is very close to the attachment probability of individual 
motors. The elasticity effects can be detected only for 
long filaments with N ^ iV*, which are softer (the effec- 
tive force constant of the filament decrease as N~^) and, 
hence, more infiuenced by the forces of the motors. For 
infinitely stiff filaments (fc — ^ oo), the crossover filament 
size diverges {N* -^ oo) and the filament elasticity is, of 
course, irrelevant on all length scales. 

Figs.[2lJA) and (B) represent two limiting cases. In the 
former, the friction is caused by the drag of the actin 
bundle in the viscous environment, while in the latter it 
originates from the attachment of the actin to the under- 



lying surface of motors. The actual friction coefficient is 
expected to lie between these two extreme values and, 
therefore, the drift velocity should exhibit an intermedi- 
ate behavior between those shown in Figs. [21^ A) and (B). 
Thus, the typical magnitude of Vdrift is expected to be of 
the order of 10 nm/scc. Interestingly, the drift velocity of 
the bundle is of the same order of magnitude as its speed 
during the bidirectional motion 10], which has also been 
found to be 2-3 order of magnitude smaller than the mov- 
ing velocity of individual myosin II motors (f ~ 6 /xm/sec 
[12|). Over a period of a few minutes the a-polar bundle 
may progress a distance of a few micrometers. This im- 
plies that the drift of a-polar bundles may be relevant to 
the active remodeling of the cell cytoskeleton occurring 
during many cellular processes. 

Our investigation of the role of the filament elastic- 
ity in modifying collective motion of molecular motors 
has been motivated by experiments which have been de- 
scribed and analyzed by using a ratchet model and a 
mean field approximation for the elastic energy [10|, llll . 
In this paper we presented a more realistic microscopic 
based model that involves the determination of the ex- 
act elastic energy of the filaments. We demonstrated 
that such a model leads to new insights and novel re- 
sults like the biased transport of filaments with no net 
polarities. Experimental verification of this surprising 
result is, however, difficult. It requires that (i) the mov- 
ing filaments are perfectly a-polar with internal (sequen- 
tial) order, and (ii) that they move for sufficiently long 
period of time such that the net drift can be extracted 
from the statistics of the unidirectional intervals of mo- 
tion. Unfortunately, the a-polar bundles are not formed 
by a well controlled process, but rather through a se- 



quence of stochastic fusion events that usually generate 
filaments with disordered, random, sequences and with 
little residual polarities ^G^. Also, in the existing experi- 
mental setup, the bidirectional motion cannot be tracked 
for more than about 10 nfinutes, which is too short for a 
meaningful statistical analysis. What should be more ex- 
perimentally testable is the other elasticity effect, namely 
the reduction in the fraction of connected motors. This 
effect, which has been attributed to the dependence of 
the elastic energy on the configuration of connected mo- 
tors (denoted by the vector C), is not limited to a-polar 
filaments. Polar filaments experiencing a non-uniform 
distribution of motor forces (i.e., when only a fraction of 
the monomers are connected to motors) will also develop 
a tensile stress that could potentially alter the attach- 
ment probability of the motors. In a future publication 
we plan to present a theoretical analysis of the attach- 
ment probability for perfectly polar filaments, similar to 
the analysis presented here for a-polar filaments. We also 
plan to investigate this effect experimentally by using a 
motility assay combined with micro-manipulation tech- 
nique (such as optical tweezers) to stall the gliding fila- 
ment and measure the mean force generated by the mo- 
tors. In the case of perfectly polar filaments, the forces of 
all the motors are applied along the same direction and, 
therefore, the total measured force should be simply pro- 
portional to the number of attached motors. 
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